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Abstract 

A search for squarks and gluinos in final states containing jets, missing transverse momentum 
and no high-p T electrons or muons is presented. The data represent the complete sample recorded 
in 201 1 by the ATLAS experiment in 7 TeV proton-proton collisions at the Large Hadron Collider, with 
a total integrated luminosity of 4.7 fbr 1 . No excess above the Standard Model background expectation 
is observed. Gluino masses below 860 GeV and squark masses below 1320 GeV are excluded at 
the 95% confidence level in simplified models containing only squarks of the first two generations, 
a gluino octet and a massless neutralino, for squark or gluino masses below 2 TeV, respectively. 
Squarks and gluinos with equal masses below 1 41 GeV are excluded. In MSUGRA/CMSSM models 
with tan p = 10, A) = and ^ > 0, squarks and gluinos of equal mass are excluded for masses below 
1360 GeV. Constraints are also placed on the parameter space of SUSY models with compressed 
spectra. These limits considerably extend the region of supersymmetric parameter space excluded 
by previous measurements with the ATLAS detector. 
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tan p = 10, Ao = and fi > 0, squarks and gluinos of equal mass are excluded for masses below 
1360 GeV. Constraints are also placed on the parameter space of SUSY models with compressed 
spectra. These limits considerably extend the region of supersymmetric parameter space excluded 
by previous measurements with the ATLAS detector. 

PACS numbers: 12.60Jv, 13.85.Rm, 14.80. Ly 



I. INTRODUCTION 

Many extensions of the Standard Model (SM) include 
heavy colored particles, some of which could be accessi- 
ble at the Large Hadron Collider (LHC) [1]. The squarks 
and gluinos of supersymmetric (SUSY) theories [2-10] 
form one class of such particles. This paper presents a 
new ATLAS search for squarks and gluinos in final states 
containing only jets and large missing transverse momen- 
tum. Interest in this final state is motivated by the 
large number of i?-parity conserving models, including 
MSUGRA/CMSSM scenarios [11-15], in which squarks, 
q, and gluinos, g, can be produced in pairs {gg, qq, qg}) 
and can generate the final state of interest through their 
direct ( q —> qXi and g — !• qqXi) and cascade decays 
to weakly interacting neutralinos, Xi, which escape the 
detector unseen. 'Squark' here refers only to the super- 
partners of the four light-flavour quarks. The analysis 
presented here is based on a study of final states which 
are reconstructed as purely hadronic. Events with recon- 
structed electrons or muons are vetoed to avoid overlap 
with a related ATLAS search [16] that requires them. 
The term 'leptons' is therefore used in this paper to refer 
only to reconstructed electrons and muons, and does not 
include r leptons. Compared to previous studies [17] , this 
updated analysis uses the full dataset (4.7 fb _1 ) recorded 
at y/s = 7 TeV in 2011 and extends the sensitivity of 
the search by selecting final state topologies with higher 
jet multiplicities. The search strategy is optimized for 
maximum discovery reach in the (rrig, mg)-plane (where 
rrig,mq are the gluino and squark masses, respectively) 
for a range of models. This includes a simplified model in 
which all other supersymmetric particles, except for the 
lightest neutralino, are given masses beyond the reach of 
the LHC. Although interpreted in terms of SUSY models, 
the main results of this analysis (the data and expected 
background event counts in the signal regions) are rele- 



vant for constraining any model of new physics that pre- 
dicts the production of jets in association with missing 
transverse momentum. 

The paper begins with a brief description of the AT- 
LAS detector (Section II), followed by an overview of 
the analysis strategy (Section III). This is followed by 
short descriptions of the data and Monte Carlo (MC) 
simulation samples used (Section IV) and of the trigger 
strategy (Section V). Section VI describes the physics 
object definitions. Section VII describes the event clean- 
ing techniques used to reject non-collision backgrounds, 
while Section VIII describes the final event selections and 
resulting event counts. Section IX describes the tech- 
niques used to estimate the SM backgrounds, with the 
systematic uncertainties summarized in Section X. Sec- 
tion XI describes the statistical model used to interpret 
the observations and presents the results in terms of con- 
straints on SUSY model parameter space. Finally Sec- 
tion XII summarizes the main results and conclusions. 



II. THE ATLAS DETECTOR 

The ATLAS detector [18] is a multipurpose particle 
physics apparatus with a forward-backward symmetric 
cylindrical geometry and nearly 4-7T coverage in solid an- 
gle [19]. The layout of the detector features four su- 
perconducting magnet systems, which comprise a thin 
solenoid surrounding inner tracking detectors and three 
large toroids used in a large muon spectrometer. Lo- 
cated between these two detector systems, the calorime- 
ters are of particular importance to this analysis. In the 
pseudorapidity region \rj\ < 3.2, high-granularity liquid- 
argon (LAr) electromagnetic (EM) sampling calorimeters 
are used. An iron/scintillator-tile calorimeter provides 
hadronic coverage over < 1.7. The end-cap and for- 
ward regions, 1.5 < |ry| < 4.9, are instrumented with LAr 
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calorimeters for both EM and hadronic measurements. 



III. ANALYSIS STRATEGY 

This analysis aims to search for the production of 
heavy SUSY particles decaying into jets and neutralinos, 
with the latter creating missing transverse momentum 
(£m? lss ). B ecause f high mass scale expected for the 
SUSY signal, the 'effective mass', m e g (defined below), is 
a powerful discriminant between the signal and most SM 
backgrounds. The requirements used to select jets and 
leptons (which are referred to as physics objects) are cho- 
sen to give sensitivity to a broad range of SUSY models. 
In order to achieve maximal reach over the (rrig^mq)- 
plane, six analysis channels are defined. Squarks typ- 
ically generate at least one jet in their decays, for in- 
stance through q — > qXi, while gluinos typically generate 
at least two jets, for instance through g — > qqXi. Pro- 
cesses contributing to qq, qg and gg final states therefore 
lead to events containing at least two, three or four jets, 
respectively. Cascade decays of heavy particles, as well as 
initial and final state radiation, tend to further increase 
the final state multiplicity. 

Inclusive analysis channels, labelled A to E and char- 
acterized by increasing minimum jet multiplicity from 
two to six, are therefore defined. In addition, the two- 
jet sample is divided into two channels, A and A', using 
the ratio of E!f lss to m e g, giving a total of six channels. 
Channel A' is designed to improve the sensitivity to mod- 
els with small supersymmetric particle ('sparticle') mass 
splittings, where the presence of initial state radiation 
jets may allow signal events to be selected irrespective of 
the visibility of the sparticle decay products. The lower 
jet multiplicity channels focus on models characterized 
by squark pair production with short decay chains, while 
those requiring high jet multiplicity are optimized for 
gluino pair production and/or long cascade decay chains. 
The final limits are set using the channel with the best 
expected sensitivity for each hypothesis. The channels 
and signal regions (SRs) are summarized in Table I. The 
final selection criteria are defined without reference to 
collision data satisfying the criteria applied earlier in the 
selection. 

The effective mass is defined to be the scalar sum of 
the transverse momenta of the leading N jets in the event 
together with £™ iss : 

N 

m cff = ]>>^+i^ iss . (1) 

i=i 

This general quantity is used to select events in two 
different ways, for which the specific values of N used 
in the sum differ. Criteria are placed on the ratio of 
E™ lss to TO ff, in which context N is defined to be the 
minimum number of jets used in the channel under con- 
sideration (for example N = 2 for channel A). In Table I, 



where the number of jets used is explicitly notated, the 
expression m e g (iVj) indicates the exact, exclusive, num- 
ber of jets used. However, the final signal selection in 
all channels uses criteria on a more inclusive definition, 
m c ff(incl.), for which the sum extends over all jets with 
Pt > 40 GeV. Requirements on m e s and E™ lss , which 
suppress the QCD multi-jet background, formed the ba- 
sis of the previous ATLAS jets + E^ iss + 0-lepton SUSY 
search [17] . The same strategy is adopted in this analysis. 

In Table I, A</>(jet,, P T miss ) mi „ is the smallest of the 
azimuthal separations between the missing momentum 
vector in the transverse plane, P T mlss , and the recon- 
structed jets. For channels A, A' and B, the selection 
requires A0(jetj, P T mlss ) m in > 0.4 radians using up to 
three leading jets. For the other channels an additional 
requirement A0(jet i; P T mlss ) min > 0.2 radians is applied 
to the remaining jets with > 40 GeV. Requirements 
on A</>(jetj, P T mlss )min and -E™ lss /m e ff are designed to re- 
duce the background from multi-jet processes. 

SM background processes contribute to the event 
counts in the signal regions. The dominant sources are: 
W+jets, Z+jets, top quark pair, single top quark, di- 
boson and multi-jet production. The majority of the 
W+jets background is composed of W — > tv events, or 
W — > ez/, [lis events in which no electron or muon candi- 
date is reconstructed. The largest part of the Z+jets 
background comes from the irreducible component in 
which Z — > vv decays generate large E™ lss . Top quark 
pair production followed by semi-leptonic decays, in par- 
ticular ti — > bbqq'rv with the r lepton decaying hadron- 
ically, as well as single top quark events, can also gener- 
ate large E™ ss and pass the jet and lepton requirements 
at a non-negligible rate. The multi-jet background in 
the signal regions is caused by poor reconstruction of jet 
energies in the calorimeters leading to apparent missing 
transverse momentum, as well as by neutrino produc- 
tion in semi-leptonic decays of heavy quarks. Extensive 
validation of the MC simulation against data has been 
performed for each of these background sources and for 
a wide variety of control regions (CRs). 

Each of the six channels is used to construct be- 
tween one and three signal regions with 'tight', 'medium' 
and/or 'loose' m e ff(incl.) selections, giving a total of 11 
SRs. In order to estimate the backgrounds in a consis- 
tent and robust fashion, five control regions are defined 
for each of the SRs, giving 55 CRs in total. Each ensem- 
ble of one SR and five CRs constitutes a different 'stream' 
of the analysis. The CR selections are optimized to main- 
tain adequate statistical weight, while minimizing as far 
as possible the systematic uncertainties arising from ex- 
trapolation to the SR, and any contamination from sig- 
nal events. This is achieved by using kinematic selections 
that are as close as possible to the relevant SR, and mak- 
ing use of other event properties to create CR samples to 
measure the backgrounds. 

The CRs are listed in Table II. CRla and CRlb are 
used to estimate the contribution of Z(—¥ i/p)+jets back- 
ground events to the SR by selecting samples of 7+jets 
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Requirement 


Channel 


A 


A' 


B 


C 


D 


E 


Trigger 
Lepton veto 


Leading jet p T > 75 GeV (EM Scale) and E^ iBB > 45 - 55 GeV 
No electron (muon) with p T > 20 (10) GeV 


££ liss [GeV] > 


160 


prO'i) [GeV] > 


130 


pr(j 2 ) [GeV] > 


60 


PtO's) [GeV] > 






60 


60 


60 


60 


P t(h) [GeV] > 








60 


60 


60 


PtOs) [GeV] > 










40 


40 


prVa) [GeV] > 












40 


A^(jeti,P T miss ) min [rad] > 


0.4 (i = {1,2,(3)}) 


0.4 (i = {1,2,3 


(■), 0.2 (pr > 40 GeV jets) 


ET se /m eS (JVj) > 


0.3 (2j) 


0-4 (2j) 


0.25 (3j) 


0.25 (4j) 


0.2 (5j) 


0.15 (6j) 


m e ff(incl.) [GeV] > 


1900/1400/- 


-/1200/- 


1900/-/- 


1500/1200/900 


1500/-/- 


1400/1200/900 



TABLE I: Criteria used to define each of the inclusive channels and streams in the analysis. The jets are ordered with the 
highest pt first. The variables used are defined in the text. The E™ 1SB /m e ff selection in any JV jet channel uses a value of m e g 
constructed from only the leading JV jets (indicated in parentheses). However, the final m e a (inch) selection, which is used to 
define the signal regions, includes all jets with pt > 40 GeV. The three m e ff(incl.) values listed in the final row denote the 
'tight', 'medium' and 'loose' selections, respectively, as used for the final SRs. 



CR 


SR Background 


CR process 


CR selection 


CRla 


Z+jets 


7+jets 


Isolated photon 


CRlb 


Z+jets 


Z(-» U)+jets 


66 GeV < m(££) < 116 GeV 


CR2 


Multi-jets 


Multi-jets 


A0(jeti, P T miss ) min < 0.2 rad 


CR3 


W(-^ fe)+jets 


W(-s> £u)+jets 


30 GeV < m T (£,Eg 1BS ) < 100 GeV, 6-veto 


CR4 


tt and single top 


tt -> bbqq'Ev 


30 GeV < m T (£,E^ iss ) < 100 GeV, 6-tag 



TABLE II: Control regions used in the analysis: the main targeted background in the SR, the process used to model the 
background, and main CR selection(s) used to select this process are given. 



and Z{— > £t)+]e%s events, respectively. The control re- 
gion CR2 uses a reversed and tightened criterion on 
A</>(jetj, P T miss )min for up to three selected leading jets 
(depending on channel) to produce a data sample en- 
riched with multi-jet background events. Otherwise it 
uses identical kinematic selections to the SRs. CR3 and 
CR4 use respectively a &-jet veto or &-jet requirement 
together with a lepton+i?™ lss transverse mass (tot) re- 
quirement to select samples of W(— > ^)+jets and semi- 
leptonic tt background events. Other selections are simi- 
lar to those used to select the corresponding signal region, 
although in CRlb, CR3 and CR4 the requirements on 
A</)(jeti, P T mlss ) m i n and E™ lss /m c g are omitted to maxi- 
mize the number of events without introducing extrapo- 
lations in energy or jet multiplicity. 

The observed numbers of events in the CRs for each 
SR are used to generate internally consistent SM back- 
ground estimates for the SR via a likelihood fit. This pro- 
cedure enables CR correlations and contamination of the 



CRs by other SM processes and/or SUSY signal events 
to be taken into account. The same fit also allows the 
statistical significance of the observation in the SR with 
respect to the SM expectation to be determined. The es- 
timated number of background events for a given process, 
iV(SR, scaled), is given by 



V(SR, scaled) = V(CR,obs) x 



N(SR, unsealed) 
iV(CR, unsealed) 



(2) 

where JV(CR, obs) is the observed number of data events 
in the CR for the process, and A^(SR, unsealed) and 
iV(CR, unsealed) are estimates of the contributions from 
the process to the SR and CR, respectively, as described 
in Section IX. The ratio appearing in the square brackets 
in Eq. (2) is defined to be the transfer factor (TF). Sim- 
ilar equations containing inter-CR TFs enable the back- 
ground estimates to be normalized coherently across all 
the CRs. The likelihood fit adjusts the predicted back- 
ground components in the CRs and SRs using the TFs 
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and the unsealed CR event counts as constraints, taking 
into account their uncertainties. The scaled values are 
output from the fit. 

The likelihood function for observing n events in one of 
the channels (A-E, loose to tight) is the product of Pois- 
son distributions, one for the signal region and one for 
each of the main control regions constraining the Z+jets 
(CRla/b), multi-jets (CR2), FF+jets (CR3) and ti (CR4) 
contributions, labelled P SR , P Z R a ,b, fjR, Avr and -Ptr 
respectively, and of the PDFs constraining the systematic 
uncertainties Cs ys t: 



L(n\fi,b, 9) 



P$K ' PzRa ■ PzRb • PjR 

■ C Syst (9). 



P^ 



WR 



Ptr 



(3) 



The total expected background is b. The expected 
means for the Poisson distributions are computed from 
the observed numbers of events in the control regions, 
using the TFs. The signal strength fi parameterizes the 
expected signal, with n = 1 giving the full signal ex- 
pected in a given model. The nuisance parameters (9) 
parameterize the systematic uncertainties, such as that 
on the integrated luminosity. 

The expected number of events in the signal region is 
denoted by As, while \ denotes the expected number of 
events in control region i. These are expressed in terms 
of the fit parameters fj, and b and an extrapolation matrix 
C (connecting background and signal regions) as follows: 



\s(n,b,0) = ,u • C S r->sr(0) • s 

+ J2 C ^M^-bj, (4) 

3 

\i{n,b,9) = pi ■ CsR^iR.(9) ■ s 

+ £<V*H(0)-&i, ( 5 ) 

i 

where the index j runs over the background control re- 
gions. The observed number of signal events in the SR 
(CRj-r) are s (bj), respectively. The diagonal elements of 
the matrix are all unity by construction. The off-diagonal 
elements are the various TFs. 

This background estimation procedure requires the de- 
termination of the central expected values of the TFs 
for each SM process, together with their associated cor- 
related and uncorrelated uncertainties, as described in 
Section IX. The multi-jet TFs are estimated using a data- 
driven technique, which applies a resolution function to 
well-measured multi-jet events in order to estimate the 
effect of mismeasurement on E>f 1BS and other variables. 
The other TF estimates use fully simulated MC samples 
validated with data (see Section IV B). Some systematic 
uncertainties, for instance those arising from the jet en- 
ergy scale (JES), or theoretical uncertainties in MC sim- 
ulation cross sections, largely cancel when calculating the 
event count ratios constituting the TFs. 

The result of the likelihood fit for each stream includes 
a set of background estimates and uncertainties for the 
SR together with a p- value giving the probability for the 



hypothesis that the observed SR event count is compat- 
ible with background alone. Conservative assumptions 
are made about the migration of SUSY signal events be- 
tween regions. When seeking an excess due to a signal 
in a particular SR, it is assumed that the signal con- 
tributes only to the SR, i.e. the SUSY TFs are all set to 
zero, giving no contribution from signal in the CRs. If 
no excess is observed, then limits are set within specific 
SUSY parameter spaces, taking into account theoretical 
and experimental uncertainties on the SUSY production 
cross section and kinematic distributions. Exclusion lim- 
its are obtained using a likelihood test. This compares 
the observed event rates in the signal regions with the 
fitted background expectation and expected signal con- 
tributions, for various signal hypotheses. Since the signal 
hypothesis for any specific model predicts the SUSY TFs, 
these exclusion limits do allow for signal contamination 
in the CRs. 



IV. DATA AND SIMULATED SAMPLES 
A. Proton-Proton Collision Data Sample 

The data used in this analysis were taken in 2011 with 
the LHC operating at a center-of-mass energy of 7 TeV. 
Over this period the peak instantaneous luminosity in- 
creased from 1.3 x 10 30 to 3.7 x 10 33 cm _2 s~ 1 and the 
peak mean number of interactions per bunch crossing in- 
creased from 2 to 12. Application of beam, detector and 
data-quality requirements resulted in a total integrated 
luminosity of 4.7 fb _1 [20, 21]. The precision of the 
luminosity measurement is 3.9%. The trigger used is de- 
scribed in Section V. 



B. Monte Carlo Samples 

MC samples are used to develop the analysis, opti- 
mize the selections, determine the transfer factors used 
to estimate the VF+jets, Z+jets and top quark pro- 
duction backgrounds, and to assess sensitivity to spe- 
cific SUSY signal models. Samples of simulated multi- 
jet events are generated with PYTHIA6 [22], using the 
MRST2007L0* modified leading-order parton distribution 
functions (PDFs) [23], for use in the data-driven back- 
ground estimates. Production of top quark pairs, in- 
cluding accompanying jets, is simulated with ALPGEN [24] 
and the CTEQ6L1 [25] PDF set, with a top quark mass 
of 172.5 GeV. Samples of W and Z/7* events with ac- 
companying jets are also produced with ALPGEN. Dibo- 
son {WW, WZ, ZZ, IU7*) production is simulated with 
SHERPA [26]. Single top quark production is simulated 
with AcerMC [27]. Fragmentation and hadronization for 
the ALPGEN samples is performed with Herwig [28, 29], 
using JIMMY [30] for the underlying event. For the 7+jct 
estimates of the Z{— > v9)+jets backgrounds, photon and 
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Z events are also both produced using SHERPA for con- 
sistency checks of the ALPGEN results. 

SUSY signal samples are generated with Herwig++ [31] 
or MadGrapri/PYTHIA6[22, 32, 33]. Signal cross sections 
are calculated to next-to-leading order in the strong 
coupling constant, including the resummation of soft 
gluon emission at next-to-leading-logarithmic accuracy 
(NLO+NLL) [34-38] [39]. The nominal cross section and 
the uncertainty are taken from an ensemble of cross sec- 
tion predictions using different PDF sets and factorisa- 
tion and renormalisation scales, as described in Ref. [40] . 

The MC samples are generated using the same param- 
eter set as Refs. [41-43] and passed through the ATLAS 
detector simulation [44] based on GEANT4 [45]. Differing 
pile-up (multiple proton-proton interactions in a given 
event) conditions as a function of the LHC instanta- 
neous luminosity are taken into account by overlaying 
simulated minimum-bias events onto the hard-scattering 
process and reweighting them according to the expected 
mean number of interactions per LHC bunch crossing. 



V. TRIGGER SELECTIONS 

The baseline triggers for the signal region event se- 
lection in the 2011 analysis use jets and E™ ms [<i6, 47]. 
The jet and E™ lss trigger required events to contain a 
leading jet with a transverse momentum (j?t), measured 
at the electromagnetic energy scale [48], above 75 GeV 
and significant missing transverse momentum. The de- 
tailed trigger specification, including the value of the 
E™ lss threshold, varied throughout the data-taking pe- 
riod, partly as a consequence of the rapidly increasing 
LHC luminosity. The trigger threshold on the missing 
transverse momentum increased from 45 GeV at the start 
of the data-taking period, to 55 GeV at the end. The trig- 
ger reached its full efficiency of > 98% for events with a 
reconstructed jet with pt exceeding 130 GeV and more 
than 160 GeV of missing transverse momentum. Trig- 
ger efficiencies are extracted using a sample selected by a 
looser trigger, taking into account correlations, i.e. cor- 
recting for the efficiency of the looser trigger. Prescaled 
single-jet triggers, which acquired fixed fractions of the 
data, are used for the trigger efficiency study. 

A second study verifies that the efficiency of the base- 
line trigger becomes maximal at the values quoted above. 
The efficiencies are determined with an independent sam- 
ple of events expected to possess E!f lss generated by neu- 
trinos. A sample triggered by electron candidates is used, 
where jets from electrons reconstructed with tight selec- 
tion criteria are discarded. This trigger selected mostly 
W — > ev events with jets and ran unprescaled, thus pro- 
viding a large number of events. 



VI. OBJECT RECONSTRUCTION 

The event reconstruction algorithms create the physics 
objects used in this analysis: electrons, muons, jets, pho- 
tons and 6-jets. Once these objects are defined, the over- 
all missing transverse momentum can be calculated. A 
failure in the calorimeter electronics created a small dead 
region (0 < r\ < 1.4, —0.8 < <fi < —0.6) in the second and 
third layers of the electromagnetic calorimeter, which af- 
fected energy measurements in about 20% of the data 
sample. Any event with a jet that is inside the affected 
region and that is expected on the basis of shower shape 
to potentially contribute significantly to the E™ ss is re- 
moved from the sample to avoid fake signals [49]. The 
energies of jets inside the affected region which are not 
expected to create E™ lss are corrected using the func- 
tioning calorimeter layers. 

Jet candidates are reconstructed using the anti-fc t jet 
clustering algorithm [50, 51] with a radius parameter of 
0.4. The inputs to this algorithm are clusters [52] of 
calorimeter cells seeded by those with energy significantly 
above the measured noise. Jet momenta are constructed 
by performing a four-vector sum over these cell clusters, 
measured at the electromagnetic scale, treating each as 
an (E,p) four- vector with zero mass. The jet energies are 
corrected for the effects of calorimeter non-compensation 
and inhomogeneities by using pt- and 77-dependent cal- 
ibration factors derived from MC simulation and vali- 
dated with extensive test-beam and collision-data studies 
[53]. Only jet candidates with p^ > 20 GeV are subse- 
quently retained. 

Electron candidates are required to have pr > 20 GeV 
and |?7| < 2.47, and to pass the 'medium' electron shower 
shape and track selection criteria described in Ref. [54]. 
Muon candidates [55, 56] are required to have matching 
tracks in the inner detector and muon spectrometer with 
p T > 10 GeV and |r?| < 2.4. 

Following the steps above, overlaps between candi- 
date jets with |7y| < 2.8 and leptons are resolved as fol- 
lows: first, any such jet candidate lying within a distance 
AR = \J (Arj) 2 + (A</>) 2 = 0.2 (cf> measured in radians) 
of an electron is discarded; then any lepton candidate 
remaining within a distance AR = 0.4 of any surviving 
jet candidate is discarded. The first requirement pre- 
vents energy deposits from being interpreted as both jets 
and electrons. The second ensures that leptons produced 
within jets are not used to veto the event during the se- 
lection described in Section VIII. 

The measurement of the missing transverse momentum 
two-vector P T mlss is based on the transverse momenta of 
all remaining jet and lepton candidates and all calorime- 
ter clusters not associated with such objects. Following 
this step, all jet candidates with | ^7 1 > 2.8 are discarded, 
owing to their lower precision. Thereafter, the remaining 
lepton and jet candidates are considered "reconstructed" , 
and the term "candidate" is dropped. 

Photons are identified with the same selection crite- 
ria as used in the ATLAS prompt photon cross section 
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analysis [57], where an isolated photon passing the tight 
photon identification criteria is required. Jets are clas- 
sified as 6-jets using a neural network algorithm, which 
takes as inputs the impact parameter measurements and 
the topological structure of &-quark decays, as described 
in Refs. [58, 59]. 



VII. REMOVAL OF NON-COLLISION 
BACKGROUNDS 

Non-collision backgrounds are produced predomi- 
nantly by noise sources in the calorimeters, cosmic ray 
events and beam collisions with residual gas in the beam- 
pipe (beam-gas events). The requirement of a vertex near 
the nominal interaction point with at least five associ- 
ated tracks is effective at suppressing these backgrounds. 
Further criteria are applied which require that the frac- 
tional energy deposited in each calorimeter layer, and 
in any cells with known quality problems, is consistent 
with that expected from beam-beam events. In addition, 
the energy observed in charged particle tracks associated 
with the calorimeter cluster, and the timing of the energy 
depositions in calorimeter cells with respect to the beam- 
crossing time are checked [53]. Following these selections, 
the remaining background is estimated by using the ob- 
served time distribution of the leading jets with respect 
to the bunch crossing, to create a background dominated 
control region. The non-collision background is found to 
be negligible in all of the SRs and CRs used. 



VIII. EVENT SELECTION 

Following the object reconstruction and event clean- 
ing described above, a lepton veto is applied to reject 
W(— > ^)+jets and leptonic tt events in which neutri- 
nos generate the £"™ ISS signature. The lepton thresh- 
old used in the veto is set at 20 (10) GeV for electrons 
(muons) to ensure that selected events correspond to a 
phase space region in which the veto efficiency is well 
understood. 

The signal regions are then defined by the kinematic se- 
lections given in Table I. Requirements on the transverse 
momenta of additional jets select inclusive 2-, 3-, 4-, 5- 
and 6-jet events in channels A/ A', B, C, D and E respec- 
tively. The jet pt thresholds for the leading up to four 
jets are set at 60 GeV in order to minimize the impact 
of pile-up on selection efficiency and improve background 
rejection. 

Removing events with a small angle in the transverse 
plane (A0) between jets and E™ lss suppresses multi- 
jet background in which mismeasurement of jet energy 
generates fake missing transverse momentum along the 
jet direction. For channels A, A' and B a requirement 
Acj) > 0.4 radians is applied to the leading (up to) three 
selected jets with px > 40 GeV, before the final SR selec- 
tion, to minimize loss of signal efficiency. For the other 



channels this requirement is augmented by a looser re- 
quirement that A(j) > 0.2 radians for all remaining se- 
lected jets with p T > 40 GeV. 

Multi-jet background is further suppressed by requir- 
ing that the E™ lss exceeds a specific fraction of the effec- 
tive mass of the event, m e ff. Coupled with the explicit 
requirement on m c g (incl.) discussed below this equates 
to a hard selection on E™ lss . The E!f ms /m, c ff value used 
decreases with increasing jet multiplicity because the typ- 
ical E^ lss of SUSY signal events is inversely correlated 
with jet multiplicity due to phase-space limitations. This 
is because additional jets in a SUSY decay chain increase 
the probability that the lightest SUSY particle (LSP) will 
be produced with low momentum through effective multi- 
body decays. Small mass splittings can also lead to low 
E™ lss . The multi-jet cross section is also suppressed at 
higher jet multiplicities, allowing the E™ 1SB requirement 
to be loosened. 

Finally, the signal regions are defined by criteria on 
TO c ff (incl.) which select events with hard kinematics in 
order to provide strong suppression of all SM background 
processes. Up to three m c s (incl.) values are specified per 
channel, corresponding to distinct signal regions 'tight', 
'medium' and 'loose', in which the final event samples 
are counted. 

Table III lists the number of data events passing each of 
the SR selections. The distributions of m e g(mcl.) (prior 
to the final m e ff(incl.) selections) for each channel for 
data and SM backgrounds are shown in Figs. 1-6. De- 
tails of the CR selections, and the methods used to obtain 
the background estimates follow in Section IX. The infor- 
mation is used in Section XI to produce the final results. 



IX. BACKGROUND ESTIMATION 
A. Introduction 

The Z(— > w)+jets process constitutes the dominant 
irreducible background in this analysis. It is estimated 
using control regions enriched in related processes with 
similar kinematics: events with isolated photons and 
jets [60] (CRla, Section IX B) and Z(-» ee/^/x)+jets 
events (CRlb, Section IXC). The reconstructed momen- 
tum of the photon or the lepton-pair system is added 
to P T mlss to obtain an estimate of the i?™ lss observed in 
Z(—t z/P)+jets events. The predictions from both control 
regions are found to be in good agreement, and both are 
used in the final fit. The small additional background 
contributions from Z{— > ee/ mjl/tt) decays in which the 
leptons are misidentified or unreconstructed, and from 
misidentified photon events, are estimated using the same 
control regions with appropriate transfer factors. The 
TF for CRla estimates Z(— > up)+}ets in the SR, and is 
corrected to give an estimate of Z+jets in the SR by mul- 
tiplying by the ratio of Z+jets events to Z{— > f P)+jets 
events derived from MC simulation. In the case of CRlb 
the TF is calculated between Z{— > ee/(U/x/Tr)+jets in 
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FIG. 1: Observed m ff(incl.) distribution for channel A. In 
the top panel, the histograms show the SM background ex- 
pectations, both before (black open histogram) and after 
(medium (red) open histogram) use of a fit to scale the ex- 
pectations to CR observations. This fit is applied to illus- 
trate the SR+CR fitting technique used in the main analysis. 
Before scaling, the MC simulation expectations are normal- 
ized to luminosity. The multi-jet background is estimated 
using the jet smearing method described in the text. Af- 
ter scaling, the VF+jets, Z+jets and ti and single top quark 
and multi-jet distributions (denoted by full histograms) are 
normalized to data in corresponding control regions over the 
full m e g range. Two MSUGRA/CMSSM benchmark model 
points with m =500 GeV, mi /2 =570 GeV, A o =0, tan ,8=10 
and fj, > and with m =2500 GeV, m 1/2 =270 GeV, A o =0, 
tan/3=10 and {i > 0, illustrating different topologies, are also 
shown. These points lie just beyond the reach of the previous 
analysis [17]. The arrows indicate the locations of the lower 
edges of the two signal regions. The bottom panel shows the 
fractional deviation of the data from the total unsealed back- 
ground estimate (black points), together with the fractional 
deviation of the total scaled background estimate from the to- 
tal unsealed background estimate (medium (red) line). The 
light (yellow) band shows the combined experimental uncer- 
tainties on the unsealed background estimates from jet energy 
scale, jet energy resolution, the effect of pile-up, the treat- 
ment of energy outside of reconstructed jets and MC simula- 
tion sample size. The medium (green) band includes also the 
total theoretical uncertainties. 

the CR and Z(—> vv/ee/ ma/ tt) -\-]ets in the SR. Thus 
both methods ultimately provide an estimate of the to- 
tal Z+jets background in the SR. 

The backgrounds from multi-jet processes are esti- 
mated using a data-driven technique based upon the con- 
volution of jets in a low E!f lss data sample with jet re- 
sponse functions derived from multi-jet dominated data 
control regions (Section IX D). Those from W+jets and 
top quark processes are derived from MC simulation (Sec- 
tion IX E). 

For each stream a likelihood fit is performed to the 
observed event counts in the five CRs, taking into ac- 
count correlations in the systematic uncertainties in the 
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FIG. 2: Observed m e fr(incl.) distribution for channel A', as 
for Fig. 1. 
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FIG. 3: Observed m e ff (fncl.) distribution for channel B, as for 
Fig. 1. 

transfer factors. 



B. .Z+jets estimate using a 7 + jets control region 

The magnitude of the irreducible background from 
Z{— > z^t>) +jets events in the SRs can be estimated us- 
ing 7 + jets data. When the vector boson pt is large, as 
required by the SR selections, the Z and 7 cross sections 
differ mainly by their coupling constants with respect to 
quarks. For this reason the cross section ratio, 

da (Z + jets)/ dp T 
Zh da (j+ jets) /dp T 1 ' 
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distribution for channel D, as 



can be used to translate the observed number of pho- 
ton events in the CR into an estimate of the number 
of Z events in the SR, taking into account the leptonic 
branching ratios of the Z and other effects. The ratio 
is expected to be robust with respect to both theoreti- 
cal uncertainties and experimental effects, related to, for 
example, jet reconstruction, which would be similar for 
both processes and therefore cancel in the ratio. 

The method uses photon events which are selected in 
two steps. The first aims to select a photon event sample 
where the efficiency and the background contamination 
are well known. The SR selections are then applied to 
these photon events, having added the photon to the 
E™ lss of the event to reproduce the E™ lss observed in 
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FIG. 6: Observed m c g (incl.) distribution for channel E, as for 
Fig. 1. 

Z(—t vv) background events. The SR selections consist 
primarily of requirements on the jets and E™ lss in the 
event, which directly or indirectly, due to the p^ recoil, 
impose kinematic constraints on the vector boson, i.e. 
the Z or photon. 

Photon events are selected by requiring at least one 
isolated photon passing the photon identification crite- 
ria discussed above. The photon trigger has an effi- 
ciency close to 100% for selected events with a photon 
Pt > 85 GeV. The photons are required to lie within the 
fiducial region \r)\ < 1.37 and 1.52 < |r?| < 2.37. After 
this first photon event selection a total of 2.8M photon 
candidates are obtained from the complete dataset, with 
an estimated purity > 95%. Figure 7(a) shows the lead- 
ing photon pt distribution for events passing the first 
photon selection. 

In the second selection step, the SR selection criteria 
from Table I are applied to the photon sample. In order 
to prevent the reconstructed photon in the event from 
also being reconstructed as a jet, jets within AR = 0.2 
of the photon are removed. The photon px is added to 
the E™ lss vectorial sum when applying the SR selections, 
using the appropriate calibration for the electromagnetic 
character of the photon shower. 

The numbers of photon candidates which are selected 
by the CRla criteria for channels A-E are presented in 
Table IV together with the numbers expected from MC 
simulation. Figure 7(b) shows the leading photon pj> 
distribution for events in CRla for SR-A medium, that 
requires m e g > 1400 GeV. Good agreement is seen be- 
tween the data and the MC simulation. 

These numbers of photons are corrected for experimen- 
tal effects as described in Ref. [57] before being used 
to estimate the TFs. The following effects are consid- 
ered. The combined identification and reconstruction ef- 
ficiency is estimated to be 86%, with an uncertainty of 
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FIG. 7: Leading photon pt distribution from data and MC simulation (a) directly after the photon selection and (b) in CRla 
for SR-A medium that requires m c g > 1400 GeV. The bottom panel shows the ratio of data to MC expectation, with the light 
(yellow) band indicating the uncertainty. 



less than 1%. The identification inefficiency is dominated 
by the tight photon identification requirements and de- 
creases with increasing photon pr- A further uncertainty 
of 5% is included to account for differences in efficiency 
of the photon isolation criteria in different event samples. 
Backgrounds from multi-jet processes and VF+jets events 
where an electron from the W decay is misidentified as a 
photon are each estimated to be ^1% for p T > 200 GeV. 
Therefore the background is neglected, but an uncer- 
tainty of 5% is assigned. 

The number of photon events selected by the CRla 
criteria is used to estimate the expected number of Z(— > 
vv) events in the corresponding SR using 



N Z ^( PT ) = iV?(pr) 

■R zh {PT)-Br{Z^vv) 



(1 - /bkg) 
£tOt) -Aripr) 



(7) 



Here _ZV 7 (pt) represents the number of photon candi- 
date events passing the CRla selections, binned in pt as 
in Fig. 7(b), /bkg the fraction of fake photons in the con- 
trol region, £ 7 (pt) the efficiency for selecting the photons 
and A y (pT) the photon acceptance. The cross section 
ratio Rz/jIpt) is determined from MC simulation. The 
uncertainties related to the cross section ratio have been 
studied using the two MC programs PYTHIA8 [61] and 
GAMB0S (an adaptation of the VECB0S program [60, 62]) 
and many of the theoretical uncertainties, such as the 
choice of scales and parton distribution functions, are 



found to cancel in the ratio, to a large extent [60]. It 
has, however, been shown that the ratio retains slight 
sensitivity to the jet selection and that multi-parton ma- 
trix elements must be used to describe correctly all the 
relevant amplitudes. The final uncertainties on R z ^(pt) 
should therefore be small, but a conservative uncertainty 
of 25% is assigned. Additional systematic uncertainties, 
common to several parts of the analysis, are discussed in 
Section X. 

The transfer factors between the CRla regions and 
their associated signal regions are obtained by averag- 
ing the correction term in the square brackets of Eq. 7 
over the measured px distribution of selected photon can- 
didates, and are given in Table V. 



C. Z+jets estimate using a Z(—>- l£) + jets control 
region 

The irreducible background from Z(— > vD)+]ets can 
also be estimated independently using the observed lep- 
tonic Z decays. The CRlb control regions are defined 
by requiring two opposite-sign electrons or muons with 
Pt > 20 GeV. In addition, the px of the leading electron 
is required to be above 25 GeV to protect against trigger 
turn-on effects. The di-lepton invariant mass must lie in 
the range 66 GeV < m(U) < 116 GeV. The E^ iss vari- 
able in the SR selection is emulated with the vectorial 
sum of the reconstructed Z boson momentum vector and 
the measured P T mlss . The SR jet and E™ ss requirements 
are applied, without selections on A(/)(jeti, P T mlss ) min or 
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SR 


Minimum Tn c a 
GeV 


7 CRla data 


7 CRla MC 
SHERPA/ALPGEN 


Est. Zuv SR (7) data 


A 


1400 
1900 


90 
9 


96 / 93.4 
9.42 / 9.33 


32.0 ± 3.4 ± 5.6 
3.2 ± 1.1 ± 0.6 


A' 


1200 


170 


176 / 180 


62 ± 5 ± 11 


B 


1900 


5 


6.21 / 6.31 


1.9 ± 0.8 ± 0.4 


C 


900 
1200 
1500 


223 
48 
6 


219 / 197 
55.8 / 44.5 
14.4 / 11.1 


64 ± 4 ± 11 
15 ± 2 ± 3 
1.9 ± 0.8 ± 0.4 


D 


1500 


3 


10.9 / 6.98 


0.86 ± 0.50 ± 0.24 


E 


900 
1200 
1400 


77 
26 
11 


71.5 / 47.4 
15.3 / 13.9 
8.71 / 6.11 


20 ± 2 ± 5 
7.7 ± 1.5 ± 1.9 
3.4 ± 1.0 ± 1.0 



TABLE IV: Numbers of photon events observed in the data and expected from the SHERPA and ALPGEN MC simulations in 
CRla for each SR, as well as the resulting estimated numbers of Z(— > vv) events in the SRs, with statistical and systematic 
uncertainties. 



E™ lss /m c ff. These changes are made to increase the ac- 
ceptance, since the precision of the method is limited by 
the rate of di-lepton events. 

In order to calculate the transfer factors, ALPGEN is 
used to estimate the number of Z+jets events in each SR 
and the number of Z{— > U)+]ets events in each corre- 
sponding CRlb control region. The uncertainties arise 
from the number of MC simulation events, the jet en- 
ergy scale and resolution, the electron and muon energy 
resolutions, the electron and muon selection efficiencies, 
the electron trigger efficiency, the electron energy scale, 
energy included in calorimeter clusters that is not associ- 
ated with physics objects, the input PDFs, the modeling 
of pile-up in the simulation, and the luminosity. 

The transfer factors themselves are listed in Table V 
and take into account the contribution from leptonic 
Z(— > tt) +jets events in CRlb. The estimated numbers 
of Z+jcts events obtained using this technique are consis- 
tent with those estimated using 7+jets events observed 
in CRla. 



D. Multi-jet background estimation 

The probability for multi-jet events to pass any of the 
SR selection cuts used in this analysis is, by design, very 
small. However, the large cross section for this process 
could potentially compensate for the low acceptance and 
hence lead to significant SR contamination. These two ef- 
fects also limit the applicability of conventional MC sim- 
ulation techniques; firstly because very large MC data 
samples are required and secondly because accurate mod- 
eling of the acceptance requires exceptionally detailed 
understanding of the performance of every component of 
the calorimeters. For this reason a data-driven method 
is used to estimate the multi-jet background in the SRs. 
The method makes use of high-statistics samples of well- 



measured data multi-jet events to minimize statistical 
uncertainties. In order to determine the acceptance of 
the SRs for poorly-measured multi-jet events, the jets in 
these events are convoluted with a function modeling the 
response of the calorimeters. This response function is 
based upon the results of MC simulations but is mod- 
ified in such a way as to give good agreement between 
multi-jet estimates and data in two additional dedicated 
analyses. This procedure minimizes the susceptibility of 
the multi-jet background estimates in the main analysis 
to systematic uncertainties arising from the Monte Carlo 
modeling of the initial response function. 

The jet response function quantifies the probability of 
fluctuation of the measured of jets and takes into ac- 
count both the effects of jet mismeasurement and con- 
tributions from neutrinos and muons in jets from heavy 
flavor decays. This function is convoluted with the four- 
vectors of jets in \ow-E™ lss multi-jet data events, generat- 
ing higher E™ ss events. These are referred to as 'pseudo- 
data' and are used to provide a minimally MC simulation 
dependent estimate of multi-jet distributions, includ- 
ing the distribution of A^Qet^, P x mlss ) m i n for high m c ff 
events. These distributions can be used to determine the 
transfer factors from the low A^Qet^, P T mlss ) min multi- 
jet control regions CR2 to the higher A<fi(jeti, P T mlss ) m in 
signal regions. 

The method, referred to as the 'jet smearing method' 
below, proceeds in four steps: 

(1) Selection of \ow-E^ lss seed events in the data. The 
jets in these events are well measured. These events 
are used in steps (3) and (4). 

(2) As a starting point the response function is 
determined in MC simulated data by compar- 
ing generator-level jet energy to reconstructed 
detector-level jet energy. 
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(3) Jets in the seed events are convoluted with the 
response function to generate pseudo-data events. 
The consistency between pseudo-data and exper- 
imental data in two analyses (see below) is then 
determined. The response function is modified and 
the convolution repeated until good agreement is 
obtained. 

(4) Jets in the seed events are convoluted with the re- 
sulting data-constrained response function to ob- 
tain a final sample of pseudo-data events. This 
sample is used to estimate the distributions of vari- 
ables defining the control and signal regions used 
in the main analysis. 

Seed events are triggered using single-jet triggers and 
offline thresholds of 50, 100, 130, 165, 200, 260 and 
335 GeV are then applied. To ensure that the events 
contain only well-measured jets, the P™ lss significance 
(defined as E% iss /y/E^ m where E T um is the scalar sum 
of the transverse energy measured in the calorimeters) is 
required to be < 0.6 GeV 1 / 2 . 

The response function is initially estimated from MC 
simulation by matching 'truth' jets reconstructed from 
generator-level particles to detector-level jets with AR < 
0.1 in multi-jet samples. The four-momenta of any 
generator-level neutrinos in the truth jet cone are added 
to the four-momentum of the truth jet. Truth jets are 
isolated from other truth jets by AR > 0.6. The re- 
sponse is the ratio of the reconstructed detector-level to 
generator-level jet transverse energy. 

A 'smeared' event is generated by multiplying each jet 
four-momentum in a seed event by a random number 
drawn from the response function. The smeared event 
E™ lss is computed using the smeared transverse momenta 
of the jets. The response function measured using MC 
simulation is modified using additional Gaussian smear- 
ing to widen the jet response, and a correction is applied 
to the low-side response tail to adjust its shape. These 
corrections improve the agreement with the data in step 

Two dedicated analyses are used to constrain the shape 
of the jet response function in step (3). The first uses 
the pt asymmetry of di-jet events. Events with two 
jets with |?7| < 2.8 and pt > 70, 50 GeV are selected, 
where there are no additional jets with \rj\ < 2.8 and 
Pt > 40 GeV. Events are vetoed if they contain any jet 
with pt > 20 GeV and n > 2.8. The pt asymmetry is 
given by 



with the -E5? 1SS in the 



Pt,i ~ PT,2 
Pt,i +PT,2 



(8) 



where the indices correspond to the jet pt ordering. This 
distribution is sensitive to the Gaussian response of the 
jets and to any non-Gaussian tails. A fit of pseudo-data 
to the collision data asymmetry distribution is used to 
adjust the response function generating the pseudo-data. 

A second analysis studies the R% distribution of > 3- 
jet events where topological selections ensure that one jet 



is unambiguously associated in 
event. The response of the detector to this jet is then 
given approximately by the quantity R2 defined by 



_ p T ■ (p{ 



P T mlss ) 



(9) 



where pj> is understood to be the reconstructed pt of 
the jet associated with the E T USS . This distribution is 
sensitive to the tails of the response function from mis- 
measured jets. When the pr of the jet is under-measured, 
P T mlss lies parallel to p T and hence R2 < 1. Conversely, 
when the pt of the jet is over-measured, P T mlss lies anti- 
parallel to p T and hence R2 > 1. Fits are performed in 
Pt and r\ bins in order to constrain the parameters de- 
scribing the low-side response function tail, which affects 
primarily the region with fl2<l. 

The i?2 distribution provides a sensitive test of the re- 
sponse function and hence of the background estimate in 
different regions of the detector, such as the transition 
between the barrel and end-cap calorimeters, where the 
energy resolution is degraded by the presence of dead ma- 
terial. The data are divided into four regions according 
to the rj of the poorly reconstructed jet associated with 
the P™ lss , shown in Fig. 8. The estimates agree well, 
with the data indicating that non-Gaussian fluctuations 
are not strongly r\ dependent. Given the good agreement 
observed between the data and estimates, no uncertainty 
is associated with the 77 dependence of the response. Fol- 
lowing this procedure, a good estimate of the jet response 
function, including non-Gaussian tails, is obtained. 

In order to illustrate the technique, Fig. 9 shows 
comparisons between SM MC simulation predictions, 
data and the jet smearing estimate for distributions 
of A0(jetj, P T mlss )min calculated with just the leading 
three jets. The figure makes use of the earlier stages 
of the event selections for SR-C loose and its associ- 
ated multi-jet control region. The final event selec- 
tions used in the analysis impose further requirements on 
A(/>(jeti,P T miss ) m i n for additional jets with p T > 40 GeV 
(see Table I). Good agreement is seen in Fig. 9 both be- 
tween the data and MC simulation, and between the data 
and the smearing estimate. 

In order to check that the above method is robust 
against changes in pile-up conditions, which changed sig- 
nificantly during data-taking, the method was repeated 
with the data divided into sub-samples corresponding 
to four time periods representative of different pile-up 
regimes. No significant dependence upon the level of pile- 
up was found. 

The resulting multi-jet transfer factors between CR2 
and SR for the signal regions are shown in Table V. 



E. W(— ► £i')+jets and tt background estimation 

The lepton veto applied to the signal events aims to 
suppress SM events with an isolated lepton. However, 
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FIG. 8: Distributions of i?2 in four bins (a-d) of |?7| of the poorly reconstructed jet, for estimated true jet px, defined as 
|pif + P T mlss | , greater than 100 GeV. The black points represent collision data while the open medium (red) histogram represents 
the combined prediction. The jet smearing method described in the text is used to estimate the multi-jet contribution (referred 
to in the plots as "pseudo-data") while MC simulation predictions are used for the other background components. The lower 
panels show the fractional deviation of the data from the prediction (black points), with the light (yellow) bands showing the 
multi-jet uncertainty combined with the MC simulation statistical uncertainty on the non-multi-jet estimate. 



such a veto does not reject all ti and VF+jets events, 
particularly when their decay products involve a lepton 
which is out of acceptance, or not reconstructed, or when 
the lepton is a hadronically-decaying r. 



To estimate the contributions from VF+jet and top 
quark backgrounds in the signal regions, two CRs are de- 
fined for each SR, one with a 6-jet veto (CR3 - enriched in 
W+jets events) and one with a &-tag requirement (CR4 
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FIG. 9: Comparison of observed and predicted distributions of A<^)(jeti, P T m 
(i = {1,2,3})), (a) after all selections except for those on A</>(jeti, P T mlss )n: 



)min for the leading three jets (A0(jeti, P T mlss ) m i n 
m e ff and E^' aa /m e s and (b) after all selections 



except for that on A<jf>(jeti, Pj 1 



for signal region C loose. The histograms show the MC simulation estimates of each 



background component. The medium (maroon) triangles show the multi-jet estimates from the jet smearing technique, normal- 
ized in the regions with A<^(jeti, P T mlss ) m in (i = {1,2,3}) < 0.2 radians, which replaces the multi-jet MC simulation estimate 
(denoted with a histogram) in the main analysis. The hatched region denotes the total uncertainty on the multi-jet estimate 
including statistical uncertainties from the seed event sample and the smearing procedure, systematic uncertainties in the jet 
response function, and bias in the seed event selection. The lower panels show the fractional deviation of the data from the 
prediction (black points), with the light (yellow) bands showing the multi-jet uncertainty combined with the MC simulation 
statistical uncertainty on the non-multi-jet estimate. 



- enriched in ti events) as defined in Table II. With the 
exception of the 6-jet requirement /veto the selections for 
CR3 and CR4 are identical and hence the two samples are 
fully anti-correlated. Both of these CRs require exactly 
one 'signal' electron or muon satisfying tighter selection 
criteria, whose transverse mass, formed with the E™ lss , 
lies between 30 GeV and 100 GeV. The lepton is then 
modeled as an additional jet, as it would be if it had 
entered the signal regions. The A<^(jet,, P T mlss ) min and 
_E™ ss /m e ff criteria which are applied in the correspond- 
ing signal regions are not applied to the CRs, in order to 
increase the CR sample sizes. 

In the electron channel, the modeling of the lepton as 
a jet is physically accurate, as the reconstruction will in- 
terpret misidentified electrons in this way. In the muon 
case, a missed muon will contribute additional missing 
transverse momentum, rather than an extra jet (although 
a small fraction of its energy may well be deposited in 
the calorimeters). When the lepton is a hadronically- 
decaying tau, the behavior lies between these two ex- 
tremes, with the hadrons being seen as jet activity and 
the r-neutrino as missing momentum. In order to be con- 
sistent between the electron and muon channels, and to 
use one high-statistics control region each for top quark 



and W events, the choice is made to model all leptons 
as jets. This is justified by the fact that the majority 
of the background comes from hadronic r-decay events, 
for which the behavior of the lepton is more jet-like than 
i?rm ss -like. n should be noted that this choice does not 
bias the background estimate because identical proce- 
dures are applied to data and to MC simulation events 
used to construct the transfer factors. The procedure has 
been validated with two alternative choices, in which the 
lepton is modeled either as missing transverse momentum 
or as a r decay. 

The transfer factors are calculated using MC simula- 
tion. Several corrections are applied to MC simulation 
events: 

• Each event in the CR is weighted by the ratio of 
the lepton identification efficiency in data to that 
in simulation. Similarly, the numbers in the signal 
region are weighted by a corresponding inefficiency 
scale factor. This weighting is performed on an 
event-by-event basis, based on the simulated lep- 
ton's transverse momentum and pseudorapidity. 

• A similar scale factor is applied for the 6-tagging 
efficiency (CR4) and fake rate (CR3), which differ 
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between data and simulation [58, 59]. This is also 
performed as an event-by-event weighting. 

• The leptons are smeared such that their energy res- 
olution reflects that measured in data. 

Various sources of systematic uncertainty on the trans- 
fer factors have been considered. For the leptons, the 
identification efficiency, energy resolution and trigger ef- 
ficiency are considered. The ^-tagging efficiency and fake 
rate, jet energy scale and jet energy resolution (for both b- 
quark and light jets separately), are considered, together 
with the effect of pile-up, of calorimeter electronics fail- 
ures and of calorimeter energy deposits not associated 
with physics objects. The fake lepton background is 
found to be negligible in both CR3 and CR4. 

The TFs between CR3, CR4, and the signal regions 
are given in Table V. Similar TFs are also computed for 
each channel between CR3, CR4 and the multi-jet control 
region CR2, where W+jets and tt events can contribute 
significantly. 

F. Estimated transfer factors 

The transfer factors estimated using the methods de- 
scribed above are summarized in Table V for each CR. 
These values, and those between the various CRs, to- 
gether with the observed event counts in each SR and 
CR form the inputs to the likelihood fit described in Sec- 
tion XI. 



X. SYSTEMATIC UNCERTAINTIES 

Systematic uncertainties arise through the use of the 
transfer factors relating observations in the control re- 
gions to background expectations in the signal regions, 
and from the modeling of the SUSY signal. For the trans- 
fer factors derived from MC simulation the primary com- 
mon sources of systematic uncertainty are the jet energy 
scale (JES) calibration, jet energy resolution (JER), MC 
modeling and statistics, and the reconstruction perfor- 
mance in the presence of pile-up. 

The JES uncertainty has been measured from the 
complete 2010 dataset using the techniques described in 
Ref. [53] and is around 4%, with a slight dependence 
upon V an d the proximity to adjacent jets. The JER 
uncertainty is estimated using the methods discussed in 
Ref. [53] . Additional contributions are added to both the 
JES and the JER uncertainties to take account of the ef- 
fect of pile-up at the relatively high luminosity delivered 
by the LHC in the 2011 run. Both in-time pile-up arising 
from multiple collisions within the same bunch-crossing, 
and out-of-time pile-up, which arises from the detector 
response to neighboring bunch crossings, are taken into 
account. 

The dominant modeling uncertainty in the MC sim- 
ulation estimate of the numbers of events in the signal 



and control regions arises from the impact of QCD jet 
radiation on m c ff. In order to assess this uncertainty, 
alternative samples were produced with reduced initial 
parton multiplicities (ALPGEN processes with 0-5 partons 
rather than 0-6 partons for W/Z+jets production, and 
0-3 instead of 0-5 for top quark pair production). 

PDF uncertainties are also taken into account. An en- 
velope of cross section predictions is defined using the 
68% confidence level (CL) ranges of the CTEQ6.6 [63] 
(including the a s uncertainty) and MSTW2008 [64] PDF 
sets, together with independent variations of the factori- 
sation and renormalisation scales by factors of two or one 
half. The nominal cross section value is taken to be the 
midpoint of the envelope and the uncertainty assigned is 
half the full width of the envelope, closely following the 
PDF4LHC recommendations [65]. 

Additional uncertainties arising from photon and lep- 
ton reconstruction efficiency, energy scale and resolution 
in CRla, CRlb, CR3 and CR4, 6-tag/veto efficiency 
(CR3 and CR4) and photon acceptance and cosmic ray 
backgrounds (CRla) are also considered. Other sources, 
including the limited number of MC simulation events as 
well as additional systematic uncertainties related to the 
response function, are included. 

Systematic uncertainties on the expected SUSY sig- 
nal are estimated through variation of the factorisation 
and renormalisation scales between half and twice their 
default values and by considering the PDF uncertain- 
ties. Uncertainties are calculated for individual produc- 
tion processes (e.g. qq, gg, etc.). 

Initial state radiation (ISR) can significantly affect the 
signal visibility for SUSY models with small mass split- 
tings. Systematic uncertainties arising from the treat- 
ment of ISR are studied by varying the assumed value 
of a s and the MadGrapri/PYTHTA6 matching parameters. 
The uncertainties are found to be negligible for large 
sparticle masses (m > 300 GeV) and mass splittings 
(Am > 300 GeV), and to rise linearly with decreasing 
mass and decreasing mass splitting to —30% for Am = 
and m > 300 GeV, and to -40% for m = 250 GeV 
and Am = 0. Signal ISR uncertainties are assumed to 
be uncorrelated with the corresponding background ISR 
uncertainties, to ensure a conservative treatment. 



XI. RESULTS, INTERPRETATION AND 
LIMITS 

The numbers of events observed in the data and the 
numbers of SM events expected to enter the signal re- 
gions, determined using the simultaneous likelihood fits 
(see Sections III and IX) to the SRs and CRs, are shown 
in Table III. The use of transfer factors between the CRs 
and SRs allows systematic uncertainties and nuisance pa- 
rameters to be dealt with in a coherent way, preserving 
any correlations, as described above. The free parameters 
are the background components in each SR, and these 
are constrained by the CR event counts and the TFs, 
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TABLE V: Summary of transfer factors from the main control regions of each background component in every stream. In CR4 
for signal regions A tight and B tight the m e ff requirements were relaxed to 1500 GeV to increase the numbers of events in the 
CRs for the minor tt background. 



within their uncertainties. The dominant irreducible 
background, from Z+jets events, is constrained by both 
CRla and CRlb, with CRla providing the largest sta- 
tistical weight. The resulting scaled predictions for the 
background components are shown in Table III. Good 
agreement is observed between the data and the SM pre- 
dictions, with no significant excesses found. The fitted 
predictions for the various background components agree 
well with the expectations from MC simulation before the 
fits, once theoretical uncertainties are accounted for. 

Data from all the channels are used to set limits on 
SUSY models, taking the SR with the best expected sen- 
sitivity at each point in parameter space. A profile log- 
likelihood ratio test in combination with the CL S pre- 
scription [67] is used to derive 95% CL exclusion regions. 
An interpretation of the results is presented in Fig. 10(a) 
as a 95% CL exclusion region in the (trig, m,j)-plane for 
a set of simplified SUSY models with m^o = 0. In these 
models the gluino mass and the masses of the squarks of 
the first two generations are set to the values shown in 
the figure, up to maximum squark and gluino masses of 
2 TeV. All other supersymmetric particles, including the 
squarks of the third generation, are decoupled. The re- 
sults are also interpreted in the tan j3 = 10, Aq = 0, /i > 
slice of MSUGRA/CMSSM models [68] in Figure 10(b). 
In these models, ISASUSY from ISAJET [69] v7.80 is used 
to calculate the decay tables, and to guarantee consistent 
electroweak symmetry breaking. 

In the simplified model with light neutralinos, with the 
assumption that the colored sparticles are directly pro- 
duced and decay directly to jets and E™ ss , the limit on 
the gluino mass is approximately 860 GeV, and that on 
the squark mass is 1320 GeV. Squarks and gluinos with 
equal masses below 1410 GeV are excluded. These val- 
ues are derived from the lower edge of the la observed 



limit band, to take account of the theoretical uncertain- 
ties on the SUSY cross sections in a conservative fash- 
ion. In the MSUGRA/CMSSM case, the limit on m 1/2 
reaches 300 GeV at high m and 640 GeV for low val- 
ues of niQ. The inclusion of signal selections sensitive 
to larger jet multiplicities has improved significantly the 
ATLAS reach at large mo- When their masses are as- 
sumed to be equal, squarks and gluinos with masses be- 
low 1360 GeV are excluded. 

In Figures 11(a) and 11(b) the limits from Fig. 10(a) 
are displayed again, but with the LSP mass set to 
195 GeV and 395 GeV respectively. For both values, 
only minor differences are seen in the limit curve, show- 
ing that the analysis retains sensitivity for a range of 
LSP masses. The signal region with the greatest reach 
is displayed at each point in the plane, showing that the 
tight, medium and loose selections all contribute to the 
final result. 

In Figure 12 limits are shown for two cases in which 
only pair production of (a) gluinos or (b) squarks is 
kinematically possible, with all other superpartners, ex- 
cept for the neutralino LSP, decoupled. This forces each 
squark or gluino to decay directly to jets and an LSP, 
as in the simplified MSSM scenario. Cross sections are 
evaluated assuming decoupled squarks or gluinos in cases 
(a) and (b), respectively. 

Similar models with only squark or gluino pair- 
production are shown in Figs. 13 and 14. However, in 
these variants, the sparticle content is augmented by an 
additional intermediate chargino with mass between the 
strongly-interacting sparticle and the LSP. This allows 
for production of additional jets or leptons and enriches 
the phenomenology. In the squark pair-production case, 
only left-handed squarks of the first and second gener- 
ations are considered in order to enhance the branch- 
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FIG. 10: The 95% CL S exclusion limits on (a) the (nig, m,})-plane in a simplified MSSM scenario with only strong production 
of gluinos and first- and second-generation squarks, with direct decays to jets and neutralinos; (b) the (mo, m^) plane of 
MSUGRA/CMSSM for tan/3 = 10, Aq — and fi > 0. Exclusion limits are obtained by using the signal region with the best 
expected sensitivity at each point. The black dashed lines show the expected limits, with the light (yellow) bands indicating 
the la excursions due to experimental uncertainties. Observed limits are indicated by medium (maroon) curves, where the 
solid contour represents the nominal limit, and the dotted lines are obtained by varying the cross section by the theoretical 
scale and PDF uncertainties. Previous results from ATLAS [17] are represented by the shaded region (blue) at bottom left in 
each case. The region excluded by chargino searches at LEP is taken from Ref. [66]. 



ing ratios of decay chains incorporating an intermediate 
chargino. The cross sections have been reduced by 50% 
to take this situation properly into account. Two dif- 
ferent parameterizations of the masses are shown. Fig- 
ures 13(a) and 14(a) vary the squark/gluino mass and the 
LSP mass, keeping the chargino mass exactly midway be- 
tween those two. In Figures 13(b) and 14(b), the LSP 
mass is instead held fixed, with the ratio of the chargino- 
LSP mass splitting to the squark/gluino-LSP mass split- 
ting defining the y-axis. When either mass splitting is 
large sensitivity to the model is enhanced by kinematics. 

The 'compressed SUSY' models suggested in Refs. [70, 
71] are also considered. In these models, the basic spar- 
ticle content and spectrum are similar to that in the 
CMSSM, but the sizes of all mass-splittings are controlled 
by a compression factor. The squark mass is set to 96% 
of the gluino mass. For presentation purposes, the lim- 
its are plotted against the gluino mass and the largest 
mass-splitting, i.e. that between gluino and LSP. Exclu- 
sion plots are shown in Fig. 15 for three classes of model: 
one in which all sparticle content is present, a second in 
which all the neutralinos and charginos apart from the 
LSP are taken to be sufficiently heavy to decouple, and 
a third in which the squarks instead are decoupled. 



XII. SUMMARY 

This paper reports a search for supersymmetry in fi- 
nal states containing high-px jets, missing transverse mo- 
mentum and no electrons with p^ > 20 GeV or muons 
with p T > 10 GeV. Data recorded by the ATLAS ex- 
periment at the LHC at y/s — 7 TeV, corresponding 
to an integrated luminosity of 4.7 fb _1 have been used. 
Good agreement is seen between the numbers of events 
observed in the signal regions and the numbers of events 
expected from SM sources. The exclusion limits placed 
on non-SM cross sections impose new constraints on sce- 
narios with novel physics. 

The results are interpreted in both a simplified model 
containing only squarks of the first two generations, a 
gluino octet and a massless neutralino, as well as in 
MSUGRA/CMSSM models with tan/3 = 10, A = 
and /x > 0. In the simplified model, gluino and squark 
masses below 860 GeV and 1320 GeV respectively are 
excluded at the 95% confidence level for squark or gluino 
masses below 2 TeV. When assuming their masses to be 
equal, squarks and gluinos with masses below 1410 GeV 
are excluded. In the MSUGRA/CMSSM case, the limit 
on toi/2 reaches 300 GeV at high mo and 640 GeV for low 
values of toq. Squarks and gluinos with equal masses be- 
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FIG. 11: The 95% CL S exclusion limits on the (ro§, m^)-plane in MSSM models with non-zero neutralino masses. Combined 
observed exclusion limits are based on the best expected CL S per grid point as for Fig. 10(a), but with an LSP mass of (a) 
195 GeV and (b) 395 GeV. Curves are as defined in Fig. 10(a). The letters overlaid on the plot show the SR that contributes 
the best sensitivity at each point. Previous results from ATLAS [17] are represented by the shaded region (blue) at bottom left 
in each case. 



low 1360 GeV are excluded in this scenario. These results 
are shown to be relatively insensitive to the assumption 
of a light LSP, up to LSP masses of about 400 GeV. 
Limits are also placed in the parameter space of a SUSY 
model with a compressed mass spectrum. 
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bands indicating the la excursions due to experimental uncertainties. Observed limits are indicated by medium (maroon) 
curves, where the solid contour represents the nominal limit, and the dotted lines are obtained by varying the cross section by 
the theoretical scale and PDF uncertainties. The letters overlaid on the plot show the SR that contributes the best sensitivity 
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